home *** CD-ROM | disk | FTP | other *** search
- #ifndef _IEEE_FP_H_
- #define _IEEE_FP_H_
-
- #include "_ansi.h"
-
- #include <machine/ieeefp.h>
-
- /* FIXME FIXME FIXME:
- Neither of __ieee_{float,double}_shape_tape seem to be used anywhere
- except in libm/test. If that is the case, please delete these from here.
- If that is not the case, please insert documentation here describing why
- they're needed. */
-
- #ifdef __IEEE_BIG_ENDIAN
-
- typedef union
- {
- double value;
- struct
- {
- unsigned int sign : 1;
- unsigned int exponent: 11;
- unsigned int fraction0:4;
- unsigned int fraction1:16;
- unsigned int fraction2:16;
- unsigned int fraction3:16;
-
- } number;
- struct
- {
- unsigned int sign : 1;
- unsigned int exponent: 11;
- unsigned int quiet:1;
- unsigned int function0:3;
- unsigned int function1:16;
- unsigned int function2:16;
- unsigned int function3:16;
- } nan;
- struct
- {
- unsigned long msw;
- unsigned long lsw;
- } parts;
- long aslong[2];
- } __ieee_double_shape_type;
-
- #endif
-
- #ifdef __IEEE_LITTLE_ENDIAN
-
- typedef union
- {
- double value;
- struct
- {
- #ifdef __SMALL_BITFIELDS
- unsigned int fraction3:16;
- unsigned int fraction2:16;
- unsigned int fraction1:16;
- unsigned int fraction0: 4;
- #else
- unsigned int fraction1:32;
- unsigned int fraction0:20;
- #endif
- unsigned int exponent :11;
- unsigned int sign : 1;
- } number;
- struct
- {
- #ifdef __SMALL_BITFIELDS
- unsigned int function3:16;
- unsigned int function2:16;
- unsigned int function1:16;
- unsigned int function0:3;
- #else
- unsigned int function1:32;
- unsigned int function0:19;
- #endif
- unsigned int quiet:1;
- unsigned int exponent: 11;
- unsigned int sign : 1;
- } nan;
- struct
- {
- unsigned long lsw;
- unsigned long msw;
- } parts;
-
- long aslong[2];
-
- } __ieee_double_shape_type;
-
- #endif
-
- #ifdef __IEEE_BIG_ENDIAN
-
- typedef union
- {
- float value;
- struct
- {
- unsigned int sign : 1;
- unsigned int exponent: 8;
- unsigned int fraction0: 7;
- unsigned int fraction1: 16;
- } number;
- struct
- {
- unsigned int sign:1;
- unsigned int exponent:8;
- unsigned int quiet:1;
- unsigned int function0:6;
- unsigned int function1:16;
- } nan;
- long p1;
-
- } __ieee_float_shape_type;
-
- #endif
-
- #ifdef __IEEE_LITTLE_ENDIAN
-
- typedef union
- {
- float value;
- struct
- {
- unsigned int fraction0: 7;
- unsigned int fraction1: 16;
- unsigned int exponent: 8;
- unsigned int sign : 1;
- } number;
- struct
- {
- unsigned int function1:16;
- unsigned int function0:6;
- unsigned int quiet:1;
- unsigned int exponent:8;
- unsigned int sign:1;
- } nan;
- long p1;
-
- } __ieee_float_shape_type;
-
- #endif
-
-
-
-
-
- /* FLOATING ROUNDING */
-
- typedef int fp_rnd;
- #define FP_RN 0 /* Round to nearest */
- #define FP_RM 1 /* Round down */
- #define FP_RP 2 /* Round up */
- #define FP_RZ 3 /* Round to zero (trunate) */
-
- fp_rnd _EXFUN(fpgetround,(void));
- fp_rnd _EXFUN(fpsetround, (fp_rnd));
-
- /* EXCEPTIONS */
-
- typedef int fp_except;
- #define FP_X_INV 0x10 /* Invalid operation */
- #define FP_X_DX 0x80 /* Divide by zero */
- #define FP_X_OFL 0x04 /* Overflow exception */
- #define FP_X_UFL 0x02 /* Underflow exception */
- #define FP_X_IMP 0x01 /* imprecise exception */
-
- fp_except _EXFUN(fpgetmask,(void));
- fp_except _EXFUN(fpsetmask,(fp_except));
- fp_except _EXFUN(fpgetsticky,(void));
- fp_except _EXFUN(fpsetsticky, (fp_except));
-
- /* INTEGER ROUNDING */
-
- typedef int fp_rdi;
- #define FP_RDI_TOZ 0 /* Round to Zero */
- #define FP_RDI_RD 1 /* Follow float mode */
-
- fp_rdi _EXFUN(fpgetroundtoi,(void));
- fp_rdi _EXFUN(fpsetroundtoi,(fp_rdi));
-
- int _EXFUN(isnan, (double));
- int _EXFUN(isinf, (double));
- int _EXFUN(finite, (double));
-
-
-
- int _EXFUN(isnanf, (float));
- int _EXFUN(isinff, (float));
- int _EXFUN(finitef, (float));
-
- #define __IEEE_DBL_EXPBIAS 1023
- #define __IEEE_FLT_EXPBIAS 127
-
- #define __IEEE_DBL_EXPLEN 11
- #define __IEEE_FLT_EXPLEN 8
-
-
- #define __IEEE_DBL_FRACLEN (64 - (__IEEE_DBL_EXPLEN + 1))
- #define __IEEE_FLT_FRACLEN (32 - (__IEEE_FLT_EXPLEN + 1))
-
- #define __IEEE_DBL_MAXPOWTWO ((double)(1L << 32 - 2) * (1L << (32-11) - 32 + 1))
- #define __IEEE_FLT_MAXPOWTWO ((float)(1L << (32-8) - 1))
-
- #define __IEEE_DBL_NAN_EXP 0x7ff
- #define __IEEE_FLT_NAN_EXP 0xff
-
-
- #define isnanf(x) (((*(long *)&(x) & 0x7f800000L)==0x7f800000L) && \
- ((*(long *)&(x) & 0x007fffffL)!=0000000000L))
-
- #define isinff(x) (((*(long *)&(x) & 0x7f800000L)==0x7f800000L) && \
- ((*(long *)&(x) & 0x007fffffL)==0000000000L))
-
- #define finitef(x) (((*(long *)&(x) & 0x7f800000L)!=0x7f800000L))
-
- #ifdef _DOUBLE_IS_32BITS
- #undef __IEEE_DBL_EXPBIAS
- #define __IEEE_DBL_EXPBIAS __IEEE_FLT_EXPBIAS
-
- #undef __IEEE_DBL_EXPLEN
- #define __IEEE_DBL_EXPLEN __IEEE_FLT_EXPLEN
-
- #undef __IEEE_DBL_FRACLEN
- #define __IEEE_DBL_FRACLEN __IEEE_FLT_FRACLEN
-
- #undef __IEEE_DBL_MAXPOWTWO
- #define __IEEE_DBL_MAXPOWTWO __IEEE_FLT_MAXPOWTWO
-
- #undef __IEEE_DBL_NAN_EXP
- #define __IEEE_DBL_NAN_EXP __IEEE_FLT_NAN_EXP
-
- #undef __ieee_double_shape_type
- #define __ieee_double_shape_type __ieee_float_shape_type
-
- #endif /* _DOUBLE_IS_32BITS */
-
- #endif /* _IEEE_FP_H_ */
-